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An experimental and theoretical study of lyotropic chromonic liquid crystals (LCLCs) confined 
in cylinders with degenerate planar boundary conditions elucidates LCLC director configurations. 

When the Frank saddle-splay modulus is more than twice the twist modulus, the ground state 
adopts an inhomogeneous escaped-twisted configuration. Analysis of the configuration yields a 
large saddle-splay modulus, which violates Ericksen inequalities but not thermodynamic stability. 

Lastly, we observe point defects between opposite-handed domains and explain a preference for point 


defects over domain walls. 

PACS numbers: 61.30.Jf, 61.30.Pq, 62.20.de, 68.37.-d 

The elastic properties of nematic liquid crystals (LCs) 
are crucial for liquid crystal display applications mm, 
and they continue to give rise to unanticipated funda¬ 
mental phenomena mu. Three of the bulk nematic LC 
deformation modes, splay, twist and bend, are well known 
and have associated elastic moduli Ki, K 2 and iFs, re¬ 
spectively. These moduli have been intensely studied be¬ 
cause they are easy to visualize, and because it is pos¬ 
sible to independently excite the modes via clever us¬ 
age of sample geometry [T 0 HT 2 ], LC boundary conditions 
[El [HI, and external fields usms]. As a result, these 
moduli have been measured for a variety of thermotropic 
and lyotropic LCs [laiisHini. By contrast, a much less 
studied fourth independent mode [2T1123] of elastic defor¬ 
mation in nematic LCs can exist; it is called saddle-splay. 
Saddle-splay is hard to visualize and to independently 
excite [ 2 SI EH. Moreover, the energy of this deforma¬ 
tion class can be integrated to the boundary, so that the 
mode does not appear in the Euler-Lagrange equations, 
and with fixed boundary conditions, the saddle-splay en¬ 
ergy will have no effect on the LC director configuration. 
Even with free boundary conditions, the saddle-splay 
energy will not affect the bulk LC configuration unless 
the principal curvatures of the surface are different, i.e., 
saddle-splay effects are not expected for spherical or flat 
surfaces. Thus, although much progress in understand¬ 
ing saddle-splay has been made [25l ES] 5 especially with 
thermotropic nematic LCs, unambiguous determination 
of saddle-splay energy effects on liquid crystal configura¬ 
tions and measurement of the saddle-splay elastic mod¬ 
ulus, 7^24, remain difficult m- 

While the bulk elastic constants described above 
strongly influence LC director configurations, LC bound¬ 
ary conditions at material interfaces also influence bulk 
structure. Indeed, considerable effort has gone into de¬ 
velopment of surface preparation techniques to produce 


particular bulk director configurations [laillHSS]. The 
saddle-splay term integrates to the boundary and effec¬ 
tively imposes boundary conditions at free surfaces favor¬ 
ing director alignment along the direction of highest sur¬ 
face curvature for positive 7^24 |34] and outwardly point¬ 
ing surface normals. For this effect to be present, the 
director cannot be held perpendicular to the surface, as 
was the case in our prior work [4]. The potential role of 
saddle-splay effects in determining bulk director configu¬ 
rations by spontaneous symmetry breaking has been ap¬ 
preciated [niiiiiT] but has been difficult to fully char¬ 
acterize; generally, molecular surface forces can impose 
preferred boundary conditions that are hard to disentan¬ 
gle from effects due to 7^24 [38l|39]. As a result, the mea¬ 
surements of 7^24 to date have wide confidence intervals 
[36ll4Ql|4T] and even vary in sign [41]. Finally, additional 
factors that have complicated assignment of saddle-splay 
effects are the so-called Ericksen inequalities [42] that re¬ 
quire 0 < 7^24 < 2 K 2 and 7^24 < 2Ki . These inequalities 
were derived assuming spatially uniform gradients of the 
director. They do not, however, apply in geometries such 
as ours in which gradients are not uniform. 

In this contribution, we investigate director configura¬ 
tions of the nematic lyotropic chromonic liquid crystal 
(LCLC) Sunset Yellow (SSY) confined within cylindrical 
glass capillaries with degenerate planar boundary condi¬ 
tions as initially reported by Refs. [43II45] . Our study 
employs a combination of polarized optical microscopy, 
measurements of director-field thermal fluctuations, and 
Frank-free-energy calculations to rationalize the observed 
structures. Importantly, we show that a large 7^24 leads 
to an escaped-twist (ET) ground state, which exhibits 
a classic double-twist configuration. Note, chiral sym¬ 
metry breaking in the ET configuration is fundamen¬ 
tally different from symmetry breaking in other LCLC 
systems with uniform principal curvatures [3], or with 
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homeotropic boundary conditions [4]. In the previous 
work, spontaneous twist deformation arises because K 2 
is much smaller than Ki and K^; K 2 A played no role in 
the energetics. In the present work, comparison of the¬ 
ory and experiment enables us to measure K 2 A for the 
first time in a LCLC; we find a value of i^ 24/^2 = 55.0, 
which strongly violates the Ericksen inequalities. Finally, 
we observe and characterize chiral hegdehog point de¬ 


fects separating chiral domains of opposite handedness. 
Interestingly, the presence of point defects rather than 
smooth domain walls also provides precise quantitative 
information about i ^24 that is consistent with our other 
conclusions. 

Before discussing the experimental results, we formu¬ 
late the theoretical problem. We assume the achiral ne¬ 
matic LCLC is described by a Frank free energy, i.e.. 


-I 


F= / d^x 


1 9 1 9 1 9 1 

(nV • n) -i^2 (n • V X n) + (n x V x n) — -i^24 V • (n x V x n -|- nV • n) 


(1) 


where n is the nematic director. Equation 0 explicitly 
includes the saddle-splay term with modulus 7^24, which 
can in principle be mimicked by a surface anchoring term 
that is coupled to surface curvature; thus we consider a 
saddle-splay term that combines the two effects [Ml [47]. 
A Rapini-Papoular type surface anchoring term with in¬ 
plane anisotropy [22[ [38] is excluded and discussed later 
in the text. The LC is contained inside a capillary of 
radius R and cylindrical coordinates are used to parame¬ 
terize its director field, n, with z along the capillary axis 
(see Fig. [^, i.e., 

n = cosa sin^d r + sinof sin/3 0 + cos^dz. (2) 

To determine the configuration of the ground state, 
we assume the director depends only on r and mini¬ 
mize the Frank free energy with respect to a{r) and 
/3(r). Degenerate planar anchoring conditions at the 
capillary surface prevent the director from having an r- 
component, so a{r = R) = 7r/2. Cylindrical symme¬ 
try sets /3(r = 0) = 0. Both a{r = 0) and /3(r = R) 
are free to vary, but stationarity of the free energy pro¬ 
vides the boundary conditions: 9^(a(r = 0) = 0 and 
RdrPir = R) = (|^ - |)sin2/3(r = R). 

With these boundary conditions, the Euler-Lagrange 
equations of the Frank free energy give [38] 


a{r) = 


TT 

2 ’ 


/3{r) = arctan 


2^K2K2iiK2i - 2K2)r/R 
VK~z[K2a - {K 24 - 2K2)ryR^]' 


(3) 

(4) 


This ET solution exists for ^^24 > 2i^2 and has right- 
handed chirality, i.e., the director streamlines form right- 
handed helices. A mirror-image solution P{r) tt —/d(r) 
exists with the same energy. Notice that the radial posi¬ 
tion r is scaled by the cylinder radius R and that Ki does 
not appear because this configuration has no splay. If 
^^24 < 2A'2 5 then only the trivial /3{r) = 0 solution exists, 
which corresponds to the simple parallel-axial configura¬ 
tion [35|. As ^^24 surpasses 2 K 2 = Kc, which is exactly 


the upper bound found by Ericksen, the system sponta¬ 
neously breaks chiral symmetry, and an ET configuration 
of one handedness grows continuously from the trivial so¬ 
lution. /3i = ;d(r = R) is plotted in Fig. Prior work 
with thermotopic LCs has found this ET configuration 
when an azimuthal anchoring condition dominates the 
behavior of I3i at the capillary surface through a chem¬ 
ical or mechanical treatment of the surface [MlEillH]. 
^{r) (Eq. 1^ can only be approximated by a linear twist 
model [34] for certain ratios of elastic constants. For LCs 
whose elastic moduli do not satisfy these ratios, such as 
SSY, polarized optical microscopy textures are strongly 
affected by the nonlinear behavior of (3{r) [47] . 

The normalized free energy of the ET configuration is 
readily calculated to be 

T = - {K 2 a - 2 K 2 ) (5) 



FIG. 1. (Color online) (a) Coordinate system used for direc¬ 
tor configuration and defect energy calculations. The director 
n is described by the angle a between the director projection 
Hxy and f, and by the angle /3 between n and the capillary 
axis, which is parallel to z. (b) 3D cutaway view of the cap¬ 
illary and the ground state director field using — 0.1, 

and K 24 IK 3 =4.6, which approximates the Frank moduli of 
30 % wt./wt. SSY at 25°C. Notice the large twist angle at 
the capillary surface close to 90°. a is independent of po¬ 
sition and /3 depends on only the radial coordinate r. This 
configuration has right-handed chirality. 
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, ^ ^K^-K2{K2^-2K2) 

H- , arctan ,_^ --, 

- K2 ^/K^{Ks + K24 - 2K2) 

where L is the length of the capillary. Notice that as K 2 A 
increases beyond 27^2, the free energy decreases contin¬ 
uously from 0, thereby confirming that the ET config¬ 
uration as a ground state is preferred over the uniform 
configuration whenever it can exist; K 2 A = 2 i ^2 marks a 
second-order phase transition line. The key to this ener¬ 
getic stabilization is the saddle-splay term: 

^ =-K2isiv? Pi- (6) 

ttL 

As noted by Ref. [34], who use an opposite surface nor¬ 
mal convention, F 24 couples the nematic director to the 
surface curvature tensor and favors alignment in the di¬ 
rection of highest curvature for K 24 > 0. In our case, this 
is the azimuthal direction along the circumference of the 
capillary. Thus, the saddle-splay free energy stabilizes 
the ET configuration despite introducing bulk director 
distortion. We also have verified that both the ET and 
the deformation-free solutions are stable whenever they 
are preferred (i ^24 > 2 i ^2 and ^"24 < 2 A' 2 , respectively); 
that is, their stability matrices have positive eigenvalues 
(see Supplemental Information) [T7] . 

Our experimental investigations used nematic SSY, a 
LCLC with relatively low twist modulus K 2 IKS ~ 0.1 
m- Briefly, five SSY samples were loaded into five dif¬ 
ferent capillary tubes with diameters 100 pm ±10%, from 
VitroCom (CV1017-100). The sealed samples were il- 



FIG. 2. Phase diagram of /3i, the angle between the ne¬ 
matic director and capillary axis at the capillary surface, as 
a function of elastic moduli ratios of saddle-splay (A 24 ) to 
twist {K 2 ) and bend (K3). Inset: an example capillary with 
streamlines indicating a surface director field at angle /3i with 
left-handed chirality. 


luminated between cross-polarizers by lOnm-bandpass- 
filtered 660 nm LED light at high (160x) magnification, 
enabling small depth of field and high spatial resolution 
imaging. Images were captured by a Uniq UP680-CL 
video camera, and a piezo-objective positioner was moved 
to image focal planes within the samples in 1 pm inter¬ 
vals. 

The capillaries without surface treatment were loaded 
with SSY and sealed to prevent evaporation. A crit¬ 
ical experimental question for any saddle-splay study 
concerns possible structure on the cylinder interfaces 
that could induce a preferred anchoring direction. To 
this end, we examined the inner capillary surfaces using 
atomic force microscopy (AEM) and scanning electron 
microscopy (SEM), and we compared the inner capillary 
surfaces to rubbed glass; the capillaries had no discernible 
grooved structures as on the rubbed glass. Since SSY is 
known to exhibit natural planar anchoring on smooth 
glass surfaces [29], our observations of the capillary sur¬ 
face strongly suggest that degenerate planar boundary 
conditions are present on the inner surfaces of the cylin¬ 
ders and any anisotropic Rapini-Papoular type anchor¬ 
ing effect would be small EH 133 113. We also considered 
alignment caused by flow during capillary filling. Load¬ 
ing capillaries with the LCLC in either the nematic or 
the isotropic phase resulted in the same type of director 
configurations. Eurther, since the filling flow is nearly 
perpendicular to the final alignment found at the cap¬ 
illary surface, flow alignment appears unlikely. Einally, 
we considered the possibility that a layer of molecules 
adsorbed to the capillary surface sets an easy access at 
the capillary surface during or shortly after filling. We 
exclude this possibility by cycling the filled capillary be¬ 
tween nematic and isotropic phase and observing that 
the director at the capillary surface retains no memory 
from cycle to cycle 07 ]. 

We measure the director angle, /d(r), directly by ob¬ 
serving a flickering speckle pattern and its direction in 
the LC. The pattern originates from director field tem¬ 
poral fluctuations and accompanying fluctuations in the 
l3i ordinary and extraordinary refractive indices which cause 
scattering [25]. These types of fluctuations of the direc¬ 
tor field have been exploited previously to measure the 
viscoelastic ratios of liquid crystals IMIl]. Our work 
follows Ref. m, which proposed using videos of LC flick¬ 
ering to discern local orientation of the director field 03. 
Elickering shape and direction depend on the local direc¬ 
tor field configuration and LC viscoelastic anisotropy. 

The experimentally measured /3{r) for one of the five 
LCLC samples studied is shown in Eig. |^. It is well fit 
by the calculated expression (Eq. 4), and the fitting pro¬ 
vides experimental values for ratios of the twist-to-bend 
and saddle-splay-to-bend elastic constants. The twist-to- 
bend ratio is in close agreement with prior measurements 
[16]. Since K 2 a/\/K^K 2 ^ tan(;di) when 7^24 > K 2 , the 
fit values become increasingly sensitive to experimental 
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FIG. 3. (Color online) Representative flickering measure¬ 
ments of 30% wt./wt. SSY in a 90.6 pm capillary at 25°C. 
(a) Single frame from a movie cropped to a 20pm square and 
after background subtraction (see main text), (b) Averaged 
FFT of a movie containing many images of fluctuations and 
contour plot of a 2D Gaussian fit to the averaged FFT. The 
long axis of the ht is perpendicular to the dominant fluctu¬ 
ation direction and yields a measurement of /3 for the im¬ 
age slice, (c) Fit of Eq. Q to /3(r) obtained by fluctuation 
measurements along the capillary radius. Error bars are the 
standard deviations in degrees of the angle found for the 2D 
gaussian hts as in (b). A nonlinear least-squares ht of the 
parameters K 2 IK 3 and K 24 IK 3 gives estimates of the elas¬ 
tic constant ratios (6.5 and 0.15, respectively) for the sample. 
Across all measurements the average K 24 IK 3 = 6.6 and has 
a bounding interval [3.8, 9.4]. 


uncertainties as /3i ^ 7r/2. For example, the ht value of 
^ 24/^3 is sensitive to the uncertainty of the measured 
size of the capillary radius, i?[47]. For the data in Fig. 1^, 
the capillary was measured to have a diameter of 90.6 pm 
to within ^ ±0.4 pm. This relatively small uncertainty, 
however, leads to the large uncertainty we give for our 
estimate of 7^24/iFs, i.e., K 24 IKS has a mean value av¬ 
eraged across experiments of 6.6 with bounding interval 
[3.8, 9.4]. By contrast, K^jK^ is a stiff parameter in the 
ht; it has a mean value averaged across experiments of 
0.12 and a standard deviation (JK 2 IK 3 = 0-04 [47] . 

We also observed hedgehog defects associated with the 
ET conhguration. In long capillaries, we typically ob¬ 
served ET domains of opposite handedness separated by 
chiral point defects. These defects were qualitatively pro¬ 
posed in Ref. [39] . We observed annihilation of neighbor¬ 
ing defects, indicating that they carry opposite topologi¬ 
cal charge [ 47 ]. The presence of nematic director singular¬ 
ities are apparent in Fig. Ef; bright-held microscopy [47] 
reveals dark spots from scattered light along the center 
of the capillary. Once found, we image the point defect 
under crossed-polarizers with the same illumination de¬ 
scribed above. We compare these experimental textures 
with those simulated numerically using Jones matrices 
m- The comparison requires a test director conhgura¬ 
tion, which we calculate using Eqs. and For conhg- 
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FIG. 4. (Golor online) (a) Left, cross-polarized quasi¬ 
mono chromatic optical images of singular point defects bor¬ 
dering ET regions of opposite handedness in 90 pm capil¬ 
lary. Yellow arrows indicate the polarizer pass axis directions. 
Right, images reconstructed using Jones matrix calculations 
from numerically computed director helds of defects, (b) A 
3D cutaway view of a capillary with opposite-handedness ET 
regions separated by a wall defect, (c) 3D cutaway view of a 
capillary with opposite-handedness ET regions separated by 
a point defect as imaged and simulated in (a). In both (b) 
and (c) the director held represents an LG with K 2 IK — 0.1 
and K24IK — 4.6, where K\ — K3 = K and the color scale is 
the same from Eig.[^. (d) Energies of the point and domain 
wall defects relative to the ET energy as a function of either 
K 24 IK or equivalently /3i, with K 2 IK = 0.1. Points indicate 
numerical calculations and lines indicate analytical approx- 
mations SZl; the latter have higher energy than the former 
but demonstrate similar qualitative behaviors. 

urations in the presence of defects, however, the direc¬ 
tor depends on both r and 2 ;; the boundary conditions 
at 2 ; ^ ±oc bring the director conhguration back to ET 
conhgurations with opposite handedness. To arrive at an 
optimized guess, we solve the Euler-Lagrange equations 
numerically with a relaxational technique 07]. The con¬ 
hgurations that emerge are very similar to what one gets 
if one takes the standard radial and hyperbolic hedgehogs 
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and simply rotates all directors by 7r/2 about the z-axis 
m- This simple operation, which is guaranteed to pre¬ 
serve hedgehog charge, automatically produces opposite 
chirality on opposite sides along ^ of the hedgehog defect 
regardless of the sign (±1) of its charge. The topologi¬ 
cal charges of successive hedgehogs necessarily alternate 
in sign [47]. Using Ki = = K, K 2 IK = 0.1 and 

K 24 /K = 4.6, numbers which are consistent with our 
measurements in the ET ground state, we observed re¬ 
markable agreement between experimental and theoreti¬ 
cal textures (Fig. [^). 

In principle, smooth domain walls can also separate 
domains of opposite handedness, in which the escaped- 
twist configuration continuously untwists from one do¬ 
main to the wall mid-plane and then continuously re¬ 
twists with opposite handedness into the other domain 
(see Fig. &) . In this case, throughout the mid-plane, 
the director would align along the capillary axis. How¬ 
ever, in SSY, we have never experimentally observed such 
a domain wall structure. Defect energetics provide an 
explanation for this observation which has an interest¬ 
ing consequence. Again, we numerically calculate the 
configurations of both domain walls and point defects to 
obtain their energies m- For these calculations, we fix 
K 2 IK = O.I in accordance with [T6| and our fluctua¬ 
tion experiments, and we allow K 24 to vary. As shown 
in Fig. [^, point defects (domain walls) have lower en¬ 
ergy than domain walls (point defects) for K 24 /K ^ 4 
{K 24 /K < 4). Using Ks = K = 6.5 pN from [16] and 
R = bO pm, a typical dimensionless energy difference of 
AF/ttRK = 0.1 corresponds to AF = 2.5 x 
where T = 298 K is the experimental temperature. If K 2 A 
is greater than the crossover value then, according 

to theory, one should not expect to observe smooth do¬ 
main walls. Thus, both our observations of defects (or 
lack thereof) and our energy analysis set 4 as an approx¬ 
imate lower bound for K 2 aIK^ in agreement with our 
fluctuation-measured value of K 2 aIK = 6.6 [3.8, 9.4]. 

In summary, we have completed an experimental and 
theoretical study of a lyotropic chromonic liquid crystal. 
Sunset Yellow, in its nematic phase and confined in a 
hollow cylinder with degenerate planar boundary condi¬ 
tions. The escaped-twist configurations found to form 
in the bulk require a large saddle-splay modulus, which 
we have measured. We also observed point defects in 
this system whose existence (compared to the absence 
of smooth domain walls) provides independent confirma¬ 
tion of i^ 24 - In the future, it will be interesting to study 
and manipulate these chiral configurations and investi¬ 
gate their formation from the isotropic phase. 
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